Qubits with electrons on liquid helium 
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We study dissipation effects for electrons on the surface of liquid helium, which may serve as 
qubits of a quantum computer. Each electron is localized in a 3D potential well formed by the 
image potential in helium and the potential from a submicron electrode submerged into helium. 
We estimate parameters of the confining potential and characterize the electron energy spectrum. 
Decay of the excited electron state is due to two-ripplon scattering and to scattering by phonons in 
helium. We identify mechanisms of coupling to phonons and estimate contributions from different 
scattering mechanisms. Even in the absence of a magnetic field we expect the decay rate to be 
< 10* s~^. We also calculate the dephasing rate, which is due primarily to ripplon scattering off an 
electron. This rate is < lO'^ s~^ for typical operation temperatures. 

PACS numbers: 03.67.Lx, 73.21.-b, 76.60.Es, 73.63.-b 



I. INTRODUCTION 

Much interest has attracted recently the idea of creat- 
ing a condensed-matter based quantum computer (QC). 
A major challenge is to have a system that vifould have 
a sufficiently long relaxation time and nevertheless could 
be controlled with high precision and allow its quantum 
state to be measured. The proposed systems include lo- 
calized electron spins in semiconductor heterostructures 
[il |, |, nuclear spins of ■^^P donors B or ^^Si nuclei 
[5| in a zero nuclear spin ^*Si matrix, electron states in 
a quantum dot excited by terahertz radiation exci- 
tons in quantum dots Jfl pi, Josephson-junction based 
systems & |l^, 0, [2[ , elect rons on helium sur- 

face |]l5| , , quantum dots coupled via a linear support 
[ p7t , and trapped polar molecules p8[ . 

The system of electrons on the surface of superfluid ''He 
is attractive, from the point of view of making a scalable 
quantum computer, because (i) it has already been ex- 
tensively studied theoretically and experimentally [^ , 
(ii) the electrons have extremely long relaxation time: 
they display the highest mobility known in a condensed- 
matter system [Q, and (iii) the inter-electron distance 
is comparatively large, ^ 1 /im. To make a QC we sug- 
gested |l^ to fabricate a system of micro-electrodes, 
which would be submerged beneath the helium surface. 
Each electrode is supposed to localize one electron above 
it, as seen in Fig. ^ and to control this electron. 

The two states of an electron qubit are the two low- 
est states of quantized motion transverse to the surface. 
To further slow down the already slow relaxation, we 
initially proposed to apply a magnetic field B± normal 
to the surface. Then the estimated relaxation time T2 
becomes as small as 10^'' s, for typical B± ~ 1.5 T and 
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FIG. 1: A sketch of a micro-electrode submerged by the 
depth h ~ 0.5 beneath the helium surface, with an elec- 
tron localized above it. The electron is driven by a field E± 
normal to the surface. This field comes from the electrode 
and the parallel-plate capacitor (only the lower plate of the 
capacitor is shown). The in-plane electron potential U^^ (r) 
is parabolic near the minimum, with curvature determined 
by the electrode potential (r = {x, y) is the in-plane position 
vector). 



temperatures T w 10 mK, whereas the clock frequency of 
the computer f2 can be in the GHz range. This attracted 
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attention of experimentalists to the project |21 
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In this paper we show that, even without a magnetic 
field, the relaxation rate of a confined electron can be 
much less than that of a free electrons. The dephas- 
ing rate can be even smaller than the previous estimate 
for a strong magnetic field. This is due to large level 
spacing in a 3D confining potential formed by a localiz- 
ing micro-electrode provided the electrode is sufficiently 
thin. Electrodes of an appropriate shape have already 
been fabricated 

For low temperatures, the major known dissipation 
mechanism is scattering by surface capillary waves, rip- 
plons. These waves are very slow. Therefore a large 
distance between electron energy levels makes it impos- 
sible to conserve energy and momentum in a one-ripplon 
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decay process. Decay of the excited electron state, i.e. 
electron energy relaxation may occur via scattering into 
two short- wavelength ripplons. We show that a very im- 
portant role is played also by decay processes, where the 
electron energy goes to phonons in helium. Such phonons 
propagate nearly normal to the surface. We identify the 
major mechanisms of electron-phonon coupling and an- 
alyze their contribution to the decay rate. Dephasing is 
due primarily to scattering of thermally excited ripplons 
off an electron. We find its temperature dependence for 
different coupling mechanisms. We also investigate the 
spectrum of sideband absorption in which a microwave- 
induced electron transition is accompanied by creation 
or annihilation of a ripplon, and analyze the related de- 
crease of the intensity of the zero-ripplon absorption line. 

In Sec. II below we analyze the energy spectrum of a 
confined electron and discuss many-electron effects. In 
Sec. Ill we discuss energy relaxation rate for different 
mechanisms of electron-ripplon and electron-phonon cou- 
pling. In Sec. IV we consider dephasing rate. Sec. V 
deals with one-ripplon sidebands and the Debye- Waller 
type factor in the zero-ripplon absorption line. In Sec. VI 
we discuss electron relaxation and dephasing from fluc- 
tuations in the underlying electrodes. Sec. VII contains 
concluding remarks. 

II. ELECTRON STATES IN ONE- AND 
MANY-ELECTRON SYSTEMS 

A. Single-electron energy spectrum 

The quantum computer considered in this paper is 
based on a set of electrons which reside in potential wells 
in free space above helium, cf. Fig. |l|. The electrons are 
prevented from penetrating into helium by a high po- 
tential barrier ~ 1 eV at the helium surface. For one 
electron, the potential well is formed by the electrostatic 
image in helium, the potential from the electrode, and 
also the potential created by the grounded plate and a 
parallel plate above the electron layer (the latter is not 
shown in Fig. |l|). 

We assume that the helium occupies the halfspace z < 
0. The image potential for an electron is — A/z, where 
A = (e-l)e^/4(£ + l), with e w 1.057 being the dielectric 
constant of helium. The energy spectrum for ID motion 
in such a potential is hydrogenic, 

En = -R/n^ in = l,2,...), R = K^m/2h^. (1) 

The effective Rydberg energy i? is w 8 K, and the effec- 
tive Bohr radius is tb = It? / Km w 76 A(to is the electron 
mass) . 

The electrode potential leads to Stark shift of the en- 
ergy levels (||) Q and to quantization of motion parallel 
to the surface. A realistic estimate of this potential and 
of the electron energy spectrum can be made by mod- 
elling the electrode as a conducting sphere with a diam- 
eter 2rei equal to the electrode diameter. The center of 
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FIG. 2: Energy difference E2 — E\ and matrix elements 
Znm ~ {n\z\m) of the electron coordinate normal to helium 
surface on the wave functions of the ground and first excited 
states of motion in the z-direction, |1) and |2), vs. the over- 
all pressing field E±. The scaled field eE^rs/R = 1 for 
E± ~ 0.91 kV/cm. The inset shows the full energy level dia- 
gram. Each level En of z-motion gives rise to a set of energy 
levels of vibrations parallel to helium surface, with typical 
spacing hu.. . 



the sphere is located at depth h beneath the helium sur- 
face. Typically we expect /i to be ~ 0.5 /zm, so that it 
largely exceeds the distance from the electron to the sur- 
face ^ Tb- For z <ti h and for the in-plane distance from 
the electrode r = {x^ + y'^f/'^ < /i, (ft.^ - r^i)^/^, the 
electron potential energy is 

;7(r,z) w-^ + ef^z + imwV^, (2) 

with 

£1. = Vc\rc\h^^ + erc\h{h? - r^j)"^, 

a;,, = {e£^/mhY/\ (3) 

Here, r — [x, y) is the electron in-plane position vector, 
and Vc\ is the electrode potential. The second term in 
£_L comes from the image of the electron in the spherical 
electrode. 

In the approximation (Q) , the electron out-of-plane and 
in-plane motions separate, with in-plane motion being 
just harmonic oscillations. Variational calculations of the 
energy spectrum of the out-of-plane motion were done 
earlier p^ . The simple model (||) with an infinite wall 
at 2; describes the observed transition frequencies 
with an error of only a few percent, which is sufficient for 
the present purposes (more realistic models have been 
discussed in literature, see Refs. |2^ and papers cited 
therein) . The full electron energy spectrum in the poten- 
tial (||) is sketched in the inset in Fig. ||. The two states of 
a qubit are the ground and first excited states of motion 
transverse to the surface, |1) and |2), both corresponding 
to the ground state of in-plane vibrations. 

In what follows, we characterize the electron state 
\i,v^mi,) with the following 3 quantum numbers: i = 
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1,2 enumerates the state of out-of-plane motion, v = 
0, 1, . . . gives the energy level of in-plane vibrations, 
and Till, = 0, 1, . . . , enumerates degenerate vibrational 
states within this IcvcL 



B. Choosing parameters of the many-electron 
system 

1. Working frequency considerations 

For a multi-qubit multi-electrode QC, the depth h by 
which the controlling electrodes are submerged into he- 
lium, the intcr-clcctrode distances , and the electrode 
potentials should be chosen in a way that would opti- 
mize performance of the QC. This includes, in the first 
place, having a high working frequency Hqc and low re- 
laxation rate T . The frequency fiqc is limited by the rate 
of single-qubit operations and by the rate of excitation 
transfer between neighboring qubits, which is determined 
by the qubit-qubit interaction. 

Single-qubit operations will be performed |l^ by ap- 
plying pulses of resonant microwave radiation, which 
cause transitions between the states |1) and |2). The 
corresponding Rabi frequency is Vl^j = e£m|zi2|/?i, where 
£m is the field amplitude. As seen from Fig. |^, \zi2\/rB ^ 
0.5, and therefore even a comparatively weak field £m = 
1 V/cm gives ^ 6 x 10^ s~^. This shows that single- 
qubit operations should not limit riqc at least at the 
level of 10^ - 10^ Hz. 

Because the wave functions of different electrons do 
not overlap, the interaction between the qubits that wc 
consider is dipolar, as in liquid-state NMR quantum com- 
puters 1^ . An important feature of electrons on helium 
is that their localization length normal to the surface 
tb greatly exceeds the atomic radius, which makes the 
dipole-dipole interaction orders of magnitude stronger 
than the dipolar interaction in atomic systems. 

Of interest to us is the part of the qubit-qubit inter- 
action that depends on the states of the qubits. Two 
types of dipole moments have to be distinguished. One 
is determined by the difference zn — 2:22 of average dis- 
tances of the electron from helium surface in the states 
|1) and |2). The dipole moment e(zii — Z22) does not 
depend on time, if we take into account time dependence 
of the wave functions, it can be called "static". The 
interaction energy between the static dipoles of the zth 
and jth qubits can be written as (l/4)?7/j*'''(T* ct^, where 
cr* = \2)i {2\i — is the operator of the difference 

of the state occupations for the ?th qubit, and 

C/J*^=e^k22-.n|V4- (4) 

The other dipole moment is associated with the 1^2 
transition. If we use time-dependent wave functions, it 
oscillates in time at high frequency ^12 = (-E2 — Ei)/h. 
Resonant interaction between such oscillating dipoles has 
energy {1/ A)U^°'"'\a\ at + H.c], where a\ = [ctL]^' = 



2|2)i (l|i is the 1 — > 2 transition operator for the ith 
qubit, and 

c4°-'=e2|zi2|V4- (5) 

The interaction between static and oscillating dipoles is 
nonresonant and can be safely neglected. 

The interactions (jj) and (||) allow implementation of a 
CNOT two-qubit gate and of interqubit excitation trans- 
fer, respectively ||l^, |l^. For a typical dipole moment 
era, the interaction energy e^r^/df^ between the qubits 
separated by dij = 1 is 2 x 10^ Hz, in frequency units. 
This energy is very sensitive to dij and can be increased 
by reducing the inter-electron distance. Eqs. (|^), (||) ap- 
ply for dij less than the distance from the electrons to the 
grounded plate in Fig. Q for larger dij the interaction is 
screened and falls down as . In practice it means that 
the interqubit coupling is likely to be limited to nearest 
and probably next nearest neighbors. 

The matrix elements Znm depend on the overall field 
E± that presses electrons against the helium surface. 
They can be obtained by solving a one-dimensional 
Schrodinger equation for the potential — Az^^ + eE±z 
with a hard wall at z = [cf. Eq. (J^; we note that the 
total field E± differs from the field £± produced by one 
electrode, see below]. The results are shown in Fig. ^. 

The difference Z22 — zn sharply decreases with increas- 
ing field for small E± because of field-induced squeezing 
of the wave functions, which is particulalry strong for 
the wave function of the excited state |2). The inter- 
play between the squeezing and better overlapping of the 
wave functions |1) and |2) with increasing field leads to 
a weak field dependence of Z12 for eE^rs/R % 1. It is 
seen from Fig. gand Eq. (Q) that, for weak pressing field 
E± < 300 V/cm, the energy of the "static" interaction is 
higher than its estimate given above by a factor varying 
from 20 to 4 with increasing E±, because of the large 
numerical value of (Z22 — zii)/rB- It is also significantly 
higher than the energy given by Eq. (^. 

In a multi-qubit system, the fields £± on different elec- 
trodes are used to tune targeted qubits in resonance with 
microwave radiation and with each other. In the simple 
case of one microwave frequency, all these fields are nearly 
the same, they differ by ^ 1 V/cm, or ~ 1%. Therefore 
in Eqs. (^), (^) we assumed that the matrix elements 
Znm are the same for different qubits. Overall, for inter- 
electron distances d ^ 1 /zm, the qubit-qubit interaction 
limits the clock frequency of the quantum computer ilqc 
to 10^ - 10^ Hz. 



2. Limitations from many-electron effects 

The electron energy spectrum should be formed so as 
to minimize the electron relaxation rate. One of the most 
"dangerous" relaxation processes is quasi-elastic scatter- 
ing by capillary waves on helium surface, ripplons, in 
which an electron makes a transition between its states 
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and a ripplon is emitted or absorbed. This scatter- 
ing is responsible for finite electron lifetime Ti. Typi- 
cal energies of appropriate ripplons are extremely small, 
~ 10"'^ K (see below). Therefore the scattering can be 
eliminated for a one-qubit system, if none of the excited 
vibrational levels of the state |1) is in resonance with the 
ground vibrational level of the state |2) shown with a 
bold line in Fig. 

From Eq. (^, for a field £± = 500 V/cm and h = 
0.5 Aim we have uj^^/2tt « 2.1 x 10^° Hz « 1.0 K. Even 
though the spacing between vibrational levels is less than 
the energy gap E2 — Ei ~ 6 — 10 K, with so big uj^^ it 
is easy to avoid resonance between E2 and an excited 
vibrational level of the state 1, i.e. between E2 — Ei and 

The situation becomes more complicated for a system 
of interacting qubits. The interaction leads to coupling 
of in-plane vibrations of different electrons. In a many- 
electron system the vibrational energy spectrum becomes 
nearly continuous. One can think that each vibrational 
level in Fig. ^ becomes a bottom of a band of in-plane 
vibrational excitations. We will assume that the width 
of the lowest band A^^ is small compared to cu^^ . The 
width of the J^th band is then ~ for not too large v. 
To avoid quasi-elastic scattering by ripplons, the electron 
energy spectrum has to be discrete, i.e. vibrational bands 
should be well separated from each other up to energies 
E2 — El, that is for v ^ {E2 ~ Ei) /Tllu^^. This means that 

A|, « hLu^^/{E2 - El). (6) 

The value of A|| depends on the geometry of the many- 
electron system. It can be found if the electrodes and the 
electrons above them form a regular 2D array, or in other 
words, the electrons form a Wigner crystal with the same 
lattice constant as the electrodes. Then, if the phonon 
frequencies of the free-standing crystal in the absence of 
the electrode potential are (k is the wave vector and 
J = 1, 2 is the branch number), then the vibrational fre- 
quencies of the pinned crystal are (tj^^ -I- oj^)^/^. The 
phonon bandwidth is small compared to lo^^ provided 
tJkj ^ in which case — maxLul^j/uj^^ ~ "^p/i^iii 

where ujp = {2Tre^nl^^ /m)^^'^ is the characteristic zone- 
boundary frequency of the free-standing Wigner crystal 
(rie is the electron density). 
It follows from the above arguments and the condition 
that quasi-elastic scattering will be eliminated for a 
pinned Wigner crystal provided 

iol^hul/{E2~Ei). (7) 

This imposes an upper limit on the nearest neighbor spac- 
ing d — mindy , because ujp cx For a square lattice 
with d = 1 fim we have LOyjlix k, 6.3 GHz. 

For the multi-electrode system, the frequency itself 
depends on the inter-electrode distance d. If the electrode 
radius rei is small compared to the depth /i, the effect 
of the electrostatic image in the electrode [in particular. 




0.4 0.7 1 1.3 1.6 1.9 

h/d 



FIG. 3: In-plane frequency Cj^^ = u^Jlo'^^ and normal to the 

surface field E\_ — E±/E'j_ for an electron above a square 
array of electrodes. The electron is localized at height h 
above one of the electrodes. The inter-electrode spacing is 
d. Electrodes are modelled by small spheres, roi//i <^ 1, 
with same positive potential The scaling frequency 

LJ|| — {eVcird/mh'^)^^^ is given by Eq. (^) and corresponds 

to the limit d — *■ 00. The scaling field is E'j_ = 2nneVcirci- 
Asymptotic behavior of uj^^ and E± for large 2nh/d is shown 
with dashed lines. 

the second term in Eq. (||) for £±] can be ignored. The 
overall potential of the electrode lattice at a distance z 
above helium surface (z -f ft, > 0) is 

V{r, z) = 2TTneVcir^i G'^ exp(iGr)e-'=(^+'*) 

-27rneV;iroi(z + h) (8) 

where G is the reciprocal lattice vector. 

The dependence of uj^^ on h/d for a square electrode ar- 
ray is shown in Fig. ^ along with the z-dependence of the 
total normal field from the electrodes. The electrostatic 
in-planc confinement is due to the spatial nonuniformity 
of the electrode potential. Therefore oj^^ falls down as 
2T:{2eVcirci/md^y^'^ ex-pi— T:h/d) for large 2'Kh/d. How- 
ever, as seen from Fig. |[ uj^^ remains close to the single- 
electrode value (||) for ft/c? < 0.5. This gives the desirable 
range of the aspect ratio h/d. 

The total perpendicular field on a localized electron E± 
comes from the electrodes and the capacitor that holds 
the system (its lower plate is shown in Fig. |^). As we 
will see, the field-induced squeezing of the electron wave 
functions (cf. Fig. ^) makes a significant effect on the 
electron relaxation rate. Therefore E± should be mini- 
mized in order to reduce relaxation effects. At the same 
time, the electrostatic confinement (the frequency ) in- 
creases with the increasing field from the electrodes. It 
would be good to compensate the out-of-plane field E± 
while keeping the in-plane potential as strongly confin- 
ing as possible. This can be accomplished using the field 
from the capacitor, which is uniform in the plane and 
does not affect in-plane confinement. 

The limitation on the compensating capacitor field 
comes from the condition that the overall field behind the 
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electron layer should attract electrons to helium, other- 
wise they will leave the surface. This field is formed not 
only by the externally applied potentials, but also by the 
electron layer itself. The total averaged over r applied 
field in the electron plane should therefore exceed Anerie. 
In other words, the uniform component 2'imeVc\rc\ of the 
electrode field —dzV\z=Q (^) can be compensated down 
to A-KeUe- The remaining pressing field on the electron 
E±_ becomes then C x 2TineVc\rc\ + Airerie with small C 
{C « 0.24 for h/d = 0.5, as seen from Fig. |). 

We note that the frequency oj^^ can be further increased 
electrostatically without increasing Ej_ by using a more 
sophisticated configuration of electrodes. Analysis of 
such configurations is outside the scope of this paper. We 
note also that, for sufficiently large lo^^ , the curvature of 
the electrode potential (||) in the z-direction may become 
substantial, particularly for highly excited states of out- 
of-plane motion. However, for a typical cu^^ /2tt = 20 GHz, 
the effective curvature-induced change of the out-of-plane 
field for lowest states 2mujfrB/e is only w 14 V/cm. 



C. Electrostatic force on helium 

Electric field from the electrodes and pressure from the 
electrons (polaronic effect) lead to deformation of the he- 
lium surface. The effect of the electrode potential can be 
easily estimated by noticing that the dielectric constant 
of helium is close to one, e — 1 sa 0.057 <C 1. Therefore 
if the surface is raised by ^(r), the associated change in 
the density (per unit area) of the free energy of helium 
AJ" in the surface field E(r) is -{e - l)El{r)^{r)/8TT. 
Bending of the surface is counteracted by surface energy, 
with density CT(9^/9r)^/2, where tr is the surface ten- 
sion. The competition between these two terms gives the 
height ^ - ^Elhd/a. For typical E^ = 3x 10^ V/cm, 
h — 0.5 /im, and d = 1 fim this gives a negligibly small 

< 10^^° cm. Therefore this effect can be safely ignored. 



III. DECAY OF THE EXCITED ELECTRON 
STATE 

A. The Hamiltonian of coupling to surface 
displacement 



density distribution. The spatial structure of surface ex- 
citations is complicated as well. However, for excitations 
with sufficiently long wavelengths to a good approxima- 
tion the vibrating helium surface can still be considered 
as a corrugated infinitely high potential wall. The elec- 
tron wave function is set equal to zero on the surface. 

In this approximation the Hamiltonian Hi of interac- 
tion of an electron with surface vibrations is obtained 
by changing the electron coordinates r ^ r, z — > z — 
^(r) where ^(r) is the surface displacement, see Refs. 
|30| , pH , p2| . The interaction is a series in the ratio £,/rB- 
Typically this ratio is very small, ^ 3 x 10~^ for thermal 
displacement with characteristic wave numbers. There- 
fore to a good approximation Hi can be expanded in ^, 



keeping only lowest-order terms. The major term, H, 
is linear in ^(r) = ^ ^qe 



(1) 



iqr 



(9) 



with 

Vpo\{x) = x^^ [1 - xKi{x 



ih 

q - — (q • P)Pz - T^q^Pz + eEi_ + Kq^Vpoiiqz), 
m Zm 

-2 



(10) 



Here, p = —ifid-r is the 2D electron momentum, and 
pz — —ihdz- The first two terms in the operator Vq de- 
scribe a kinematic interaction, which arises because the 
electron wave function is set equal to zero on a non-flat 
surface. The term Vpoi{qz) describes the ch ang e of the po- 
larization energy due to surface curvature |^lf |32| {Ki{x) 
is the modified Bessel function). 
The quadratic in ^ coupling is 

^i^^ = Yl ^qi'^q2 cxp[i(qi + q2)r] V^q^q^. (11) 
"^^qi ,q2 

As in the case of linear coupling, it also has kinematic 
P2[ and polarization parts. 



q.q. = + VqM. = -(qiq2)p^/2m, (12) 



with 

V;(P°i) - -Az-3[1 - uiqiz) - u{q2z) + u(|qi + qajz)] 



u{x) = x^K2{x)/2 



(13) 



The major mechanism of electron relaxation for low 
temperatures is scattering by vibrations of the liquid he- 
lium surface. A complete calculation of the energy of 
coupling to surface vibrations is nontrivial. The density 
profile of the interface between helium and its vapor has 
a complicated form, with the 10%/90% interfacial width 
« 6 — 7 A for low temperatures p9| . As a consequence, 
even for a flat surface the electron potential is more com- 
plicated than the simple image potential —A/z for z > 
and a sharp wall at z = (^ |2^. In particular the re- 
pulsive barrier is smooth, but it becomes high compared 
to the binding energy R already on the tail of the helium 



1. Coupling to ripplons 

The biggest contribution to surface vibrations comes 
from capillary waves, ripplons. The displacement ^q is 
related to the creation and annihilation operators of rip- 
plons by 

^q = S-'^Hhq/2pu;,)^^^{bq + bl^), 

where S is the area of the system, p is the helium den- 
sity, and the ripplon frequency cjq = {aq'^ / p)^/'^ for 
q > {pg/aY/"^. 
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The change of variables used to take into account the 
hard waU potential on helium surface leads also to extra 
terms in the kinetic energy of ripplons coupled to the elec- 
tron, which is yet another source of electron-ripplon cou- 
pling Compared to similar terms in Eqs. (^), (11), 
these terms have an extra parameter LUqin/hq^, which is 
extremely small for typical q. 

There are several limitations on the wave numbers q 
of ripplons for which the electron-ripplon coupling has 
the form (|), (|ll|). Monarkha and Shikin ||| argue that 
essentially qrs should be ^ 1. Clearly, q should be small 
compared to the reciprocal width of the helium liquid- 
vapor interface and the reciprocal decay length of the 
electron wave function into helium (note that there is no 
factor 27r here, because a capillary wave with wave num- 
ber q decays into helium as exjp{qz), for a sharp inter- 
face) . Both lengths are of order of a few angstroms, which 
means that the large-g cutoff (/max should be ^ 10^ cm^^. 

A cutoff at 10^ cm~^ is consistent also with the condi- 
tion that H^'^'' ( pi] ) be small. To first order in the kine- 

(2) 

matic part of , which dominates for large q, the rel- 
ative change of the electron kinetic energy for motion 
transverse to the surface for T = is 

SK/K = hqUl/Uniap)'/'. 

This gives 6K/K w 3 x 10"" for q^ax = 10^ cm^^ (for 
<Zmax = 10^ cm"^ the correction would be equal to 1). 
Presumably the interaction with ripplons as a whole and 
in particular inelastic scattering by ripplons falls down 
for q much bigger than 10^ cm~^, because for such mo- 
mentum transfer an electron "resolves" atomic structure 
of helium. Finding Hi for such q requires a full calcu- 
lation of the ripplon-induced modulation of the electron 
potential for the diffuse helium surface, which is not the 
subject of the present paper. 

In what follows we will use spectroscopic notations and 
evaluate the decay rate as F = l/2Ti, where Ti is the life- 
time of the excited state. Defined in this way, F gives the 
decay rate of the off-diagonal matrix element pi2 of the 
electron density matrix and the decay-induced broaden- 
ing of the absorption line. 



B. One-ripplon decay 

The important consequence of strong in-plane confine- 
ment is that it essentially eliminates decay processes in 
which an electron transition is accompanied by emission 
or absorption of a ripplon. This happens because ripplons 
are very slow, and energy conservation in a transition re- 
quires transfer of too large momentum for an electron to 
accommodate. 

For low temperatures, fc^T <C hto^p qubit relaxation 
is due to electron transitions |2,0, 0) — > \l,iy,m^) from 
the ground vibrational level of the state |2) of z-motion 
into excited states of in-plane vibrations in the state |1) 
of z-motion, see Fig. 0. Minimal energy transfer is of the 



order of the vibrational level spacing hu^^ . It corresponds 
to a transition into vibrational states with the energy 
closest to E2 from below, i.e. with v = I'c ^ int[(£^2 ^ 
Ei)/huj^^]. The squared matrix element of the transition 
Kv = rriu = 0| exp(iqr)|z/c,TO,y^)p is exponentially small 
for q ^ i^y^/aii , where 



= (aW)V2 



(14) 



is the electron in-plane localization length. 

1/2 

The frequency of ripplons with q — Vc /ay is much 
less than uj^^ provided 

c.,, :^{a/pY'^[m{E2-E^)/h^]^'\ 

This inequality is satisfied already for lo^J^tt > 0.2 GHz, 
whereas a typical lo^^ for confined electrons is ^ 20 GHz. 
Therefore one-ripplon decay is exponentially improbable. 
This result does not change for a many-electron system 
provided the bands of in-plane vibrations are narrow, as 
discussed in Sec. II B, see Eq. (||). 



C. Two-ripplon decay 

Even for large separation between electron energy lev- 
els, where one-ripplon decay processes are exponentially 
suppressed by the restriction on the transferred mo- 
mentum, decay into two ripplons may still be possible 
p^ , Hj. Indeed, each of the wave vectors qi,q2 of the 
participating ripplons can be big, this is only their sum 
that is limited by the typical reciprocal electron wave- 
length. The energy and momentum conservation law re- 
quire than that the ripplons have nearly same frequen- 
cies and propagate in opposite directions, qi « — q2 and 



^91 



"92 



5E/2h, where 5E is the change of the elec- 
tron energy. 

A typical minimal value of 5E for the decay of the elec- 
tron state 1 2) is determined by the distance between the 
energy levels of in-plane vibrations fiuj^^ . The ripplon fre- 
quency LOq/2-K becomes equal to typical uj^^ /At: = 10 GHz 
for q = 1.2 X 10^ cm~^, i.e. presumably beyond the 
range of applicability of the electron-ripplon coupling 
theory (|), (0). Therefore we use this theory to esti- 
mate the rate of decay only with minimal energy trans- 
fer, i.e. we are interested in ripplon-induced transitions 
|2, 0, 0) ^ |1, mn) with v = Vc- In this case 

5E = E2 — El — Vchio^^ ■ 

In calculating the total rate of decay into different 
states of the vih vibrational energy level of a two- 
dimensional oscillator it is convenient to use the relation 

= x''e^yiyl, x^q^a^^/2 (15) 
independent of j — 1,2. 
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FIG. 4: Scaling factors K12 ® and Ui2{q) (0) in the es- 
timated probabilities of scattering into two ripplons due to 
the kinematic and polarization couplings, respectively. The 
ripplons propagate in opposite directions with nearly same 
wave numbers ijrcs given by the energy conservation condi- 
tion 2ijjq = SE/h for q — Qics- Only transitions with smallest 
energy transfer 5E are taken into account and the approxi- 
mation of an infinite sharp potential wall for an electron on 
helium surface is used. The data for U12 refer to gros ~ 3.5/rs , 
which corresponds to 5E/2-kTi ~ 5 GHz. 



For the kinematic two-ripplon coupling, which is given 
by Eq. (0) with Km = 
account taken of (|lj) is 



^qfqa : the dccay rate Fg^^ with 



r(k) _ 



2 7/2 



247ra2pi/2a3/2 ' 



Ki2 = 



(l|P^/2m|2) 
R 



(16) 



place K2{\qi + q2|z) with its small-z limit in (jlj). Then 
Ui2{q) - 2rl{l\z-^[2 ~ q^z^K2{qz)]\2). (17) 



The coefficient U12 as given by Eq. (17) is shown in 



Fig. 1^. For 5E and lo^^ chosen above we haveF^r'Vr;^'^ 
0.1 for £;_L = - 300 V/cm. The rate F^p°'^ grows much 

slower then r2^.^ with increasing g^cs (and thus with in- 
creasing 6E). 



Besides the terms Fj^^ and T)^"'' due to purely kine- 
matic and polarization mechanisms, there is a contribu- 
tion to the decay rate from the interference of these two 



(pol) 



mechanisms. It is smaller then F 
be discussed. 



(k) 



F 



(pol) 



and will not 



D. Phonon-induced decay 

An important channel of inelastic electron scattering is 
decay into phonons in helium. For a typical energy trans- 
fer 5E huj^^ , the wave numbers of the phonons partici- 
pating in the decay are 'j-'h /^s {vs is the sound velocity 
in helium). On the other hand, the in-plane momentum 
transfer is limited to ~ ^/ay ^ Tt-ld^Jvs- As a result, 
only phonons propagating nearly normal to the surface 
(in the z-direction) may be excited in a one-phonon decay 
(cf. Ref. g). 



Decay due to phonon-induced surface displacement 



Here, ^rcs is the ripplon wave vector given by the energy 
conservation law, ojq^^^ — 5E/2h, and we assumed that 

1 /2 

Qics ^ Vc /«![, which is the condition for the ripplons 
created in the transition to propagate in opposite direc- 
tions. The scaled matrix element of the kinetic energy 
if 12 is shown in Fig. IJ. 

The rate T^2r depends on ^los and therefore on 5E 
very steeply, T^^r 5E'^I^. For 5E = fu^^J^ and 
io^j2'K = 20 GHz we have gj-cs ~ 4.6 x 10^ cm^^ and 

F^^^ = 7.6 X 10^ - 3.8 X 10^ s^^ for the pressing field 
E± = — 300 V/cm. This value can be decreased by 
reducing 5E. However, Eq. (^6|) is probably an overes- 
timate even for the 5E used above, because it is based 
on the approximation of an infinite-wall potential for an 
electron and the assumptions that the helium surface is 
sharp. 

The expression for the two-ripplon decay rate F2^.°''' 
due to the polarization two-ripplon interaction (|l3| ) has 

the same form as Eq. ([ij). Just the factor if^g^r/s^ in 

Eq. ( p!6| ) has to be replaced with Ui2iqi:es)/r%qrcs , where 

U12 is determined by the matrix element of 2r'^Vqfq2 
(|^) on the functions |1), |2). The major contribution to 
this matrix element comes from the range of compara- 
tively small z, and therefore for an estimate one can re- 



Two coupling mechanisms are important for decay into 
phonons. One is related to phonon-induced displace- 
ment of the helium surface. This mechanism can be 
quantitatively described in the approximation of a sharp 
helium boundary, which provides an infinitely high po- 
tential barrier for electrons. The coupling is given by 
Eqs. (^), ( p^ ) with being a phonon-induced compo- 
nent of the surface displacement. As in the case of cou- 
pling to ripplons, it would be unreasonable to use it for 
short-wavelength phonons, in particular for phonons with 
q, > 10^ cm^i. For typical u;J2tt = 20 GHz we have 
Qz ^ /vs ~ 5 X 10^ cm~^. Therefore we will again con- 
sider decay of the state |2,0, 0) into closest lower-energy 
states |1, J^, m^). 

For typical q = {qx.qy) ^ l/ay and < u^Jvs we 
have aq^ / pv^qz <C 1. This inequality allows one to 
think of helium surface as a free boundary for phonons 
and to ignore coupling between phonons and ripplons 
p6| , ISTj . Then surface displacement is simply related to 
the Fourier components uq of the phonon displacement 
field [here, Q — (q, qz) is the 3D phonon wave vector, and 
uq is the displacement along Q]. In turn, uq is related 
to the operators of creation and annihilation of phonons 
in a standard way. 



{n/2pVvsQY'\ 



(18) 
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{V is the volume of helium) . 

With @, we obtain the rate of decay |2,0, 0) 
\l,Vc,mi,) due to phonon- induced surface displacement 

-'^ph form 



m„=0 



dq|(2,0,0|e^'i%|l,j/„m,)|2.(19) 



is] 

We start with the contribution to Fpj^ from the kine- 
matic terms in [the first two terms in Eq. (|l^)]. Tak- 
ing into account that the diagonal with respect to out- 
of-planc motion matrix element 

(j,0,0|exp(zqr)[(qp) + V/2]b- 

= -fmuj^^ {j, 0, 0| exp(iqr)|j, ly, m^) {j = 1, 2) 

and using Eq. (|l5|), we obtain for the kinematic contri- 
bution 



ph 



iE2 



12 



Airpvsh'^SE 



(20) 



The numerical value of T 



(s;k) 
ph 



is 7.8 X 10^ s-i for Si = 0, 



io^^jl-K ~ 21.1 GHz, and ^E k. huj^^ {vc = 5 in this case). 
It goes up to - 1.5 X 10^ s-i for Ei_ = 300 V/cm and 
a;||/27r = 20.6 GHz (in this case Vc — 12). The values 
of were adjusted here to meet the condition 5E — 
E2—Ei — Vifi^\\ ~ ff-^ii for the energy spectrum calculated 
for a sharp helium boundary; the real level spacing is 
a few percent smaller |2j, |2^, |2^, leading to a slightly 
smaller v'^^^^ for ujj2Tr ~ 20 GHz. We expect a more 



significant change (reduction) of T^p^^'' due to diffuseness 
of helium surface. 

(s) 

The contribution to Epj^ from the polarization term in 
01)1 has the form 



Vq [the last term in Eq, 



n(s:pol) 

ph 



4i?2r| 



vj T^pvs SEaf Jq 



dxe~''x''''+^v\x), (21) 



where v'^{x) = |(l|wpoi[(2x)-^/^2;/a||]|2)|^. The numerical 
value of r^''j^''°'^ is ~ 7 X 10^ s~^ for E_i_ =Q and goes up 
to ~ 7 X 10^ s"^ for E±_ — 300 V/cm (we used same uj^^ 
as in the above estimate of ^'^^^f' ) . 

There exists also a contribution to T^^^ ( p^ ) from 
the interference term, which is bilinear in the polariza- 
tion and kinematic interaction energies ( p^ ) (the terms 
■'^ph^'^"''' quadratic in these interactions). It can be 
obtained from (|9|) in the same way as Tp^j^'^'^"'''. It is 
positive and of the same order of magnitude as Pp'^jj'^'''"'''. 

The scattering rate Epj^ can be reduced by decreasing 
the pressing field E^. More importantly, it can also be 



reduced by increasing the frequency uj^^ , in which case the 
wavelength of the phonons to which the energy is trans- 
ferred will become smaller then the width of the diffusive 
layer on helium surface (in fact, the above calculation 
probably already overestimates the scattering rate). A 
transition to such frequency can be accomplished with a 
magnetic field applied transverse to the surface, as ini- 
tially suggested in Ref. O. 



2. Decay due to phonon-induced modulation of the helium 
dielectric constant 

Another mechanism of coupling to phonons is through 
phonon-induced modulation of the image potential of an 
electron. It results from the modulation of the helium 
density 5p and related modulation of the dielectric con- 
stant 5e. It is reasonable to assume that, for long wave- 
length phonons, 5e = {e — \)5p/ p. To lowest order in 
£ — 1, (5e the coupling energy is 



(d) 



-!- / dR'fe(R')£;2(R';R) 
Stt j 



Here the integration goes over the space occupied by he- 
lium, R = (r, z) is the 3D position vector, and E(R'; R) 
is the electric field at R' created by an electron located 
at a point R. 

The coupling Hamiltonian can be written in the form 



UQexp(zqr)V(^''\ V^'^ = iKqQv'^''^ (22) 



(d) 



with u^'^) = v'''^\q,qz, z) being 



„(d) 



dz'{z + z')"^e-'«-^'ii:i [q{z + z')] . (23) 



As in the case discussed in the previous section, the 
coupling (|2^) gives rise to phonon-induced electron tran- 
sitions between electron energy levels accompanied by 
emission of a phonon. Here, too, typical in-plane wave 
numbers of emitted phonons q are much less than the 
normal to the surface wave number « SE/hu. The 
expression for the corresponding decay rate has the 
form 



-(d) 
ph 



SErl 
nh pVg 



rfgg3|(2|^,(d)|i)|2g(^^,5) (24) 



Evaluation of the integral is largely simplified by the fact 
that the function q^g{i', q) sharply peaks at g = « 
(2j/ -I- 3)i/Va||- Therefore, with an error less than 10% 
one can replace w''^' in (j2^) by its value (p3|) for q — q^^ 

For u;||/27r = 20GHz and 6E = huj^^ 
varies from ^ 1 x 10'' s^^ to ~ 6 x 10"* s" " wim £j± 
increasing from to 300 V/cm. However, these values 
have to be taken with care. The integrand in v^'^^ (23) is 



the value of Fp"^^ 
"1 with E 



a fast oscillating function of z' on the characteristic scale 



9 



z' ^ tb, because typically qzrs ^ 1 (^z^b « 4 for cho- 
sen In addition, the matrix element of w*-'^-' in ( |2^ ) 
has an integrable singularity for z = z' = (the wave 
functions ipniz) oc z for z 0). As a result, a significant 
contribution to the matrix element comes from small dis- 
tances from the helium surface, z' ^ tb- Changing, in 
view of diffuseness of helium surface, the limit of integra- 
tion in (|2^ ) from z' = to a more reasonable z' — rs/lO 
reduces the value of rp"^^"* by a factor of 3. 

The decay rate v'"^ decreases with the increasing SE 
roughly as 1/5E (and even faster, in view of the "dead" 
layer on the diffuse surface). For higher SE and, re- 
spectively, for higher wave numbers of resonant phonons, 
the simple approximation ( p2[ ) no longer describes the 
electron-phonon interaction. Therefore, as in the case of 
scattering due to phonon-induced surface deformation, a 
way to reduce the scattering rate is to increase the fre- 
quency of in-plane vibrations. 

Full coupling to phonons is given by the sum of all 
couplings discussed in this section, with Vq in Eq. ( |l9| ) 

replaced by Vq -I- aciVQ^' with aq w 1 for typical Q. 
The total rate of scattering by phonons contains cross- 
terms which describe interference of different coupling 
mechanisms. As mentioned above, we omit these terms, 
because they do not change the overall estimate of the 
rate. We note that an interesting situation may occur if 
one of the transition frequencies of the electron comes in 
resonance with the roton energy. In this case we expect 
an increase of the decay rate. Observing it would be a 
direct demonstration of coupling to volume excitations 
in helium. 



IV. DEPHASING DUE TO RIPPLON 
SCATTERING 

In addition to depopulation of the excited state of a 
qubit, electron coupling to excitations in liquid helium 
leads also to decoherence or dephasing, i.e. to decay of 



the phase difference between the qubit states |2, 0, 0) and 
|1, 0, 0). The mechanism of this decay is random modu- 
lation, by thermal fluctuations in helium, of the distance 
between the energy levels 1 and 2. In other terms it 
can be described as quasi-elastic scattering of thermal 
excitations off an electron. The scattering is different 
in different electron states. Therefore it randomizes the 
phase difference between the wave functions of the states 
without causing interstate transitions. The correspond- 
ing mechanism is known for defects in solids [3^ as mod- 
ulational or Raman broadening. For electrons on helium 
it was discussed in Refs. |l5|, 

Dephasing comes primarily from coupling to ripplons, 
because these excitations are soft. A typical wave number 
Qr and frequency cOr of ripplons coupled to an electron are 
given by 

qr = , LOr = UJ,^ = [a/pY'\l'\ (25) 
For ujJ2ti = 20 GHz we have ujr/2n « 4.8 x lO"^ Hz 
« 2.3 mK. Therefore even for temperatures as low as 
10 mK ripplon occupation numbers are large. To the 
lowest order of perturbation theory, quasi-elastic ripplon 
scattering by an electron is determined by two-ripplon 
coupling, with the Hamiltonian 

i=l,2q,q' 

Individual terms in the sum over q, q' describe scattering 
of a ripplon with wave number q' into a ripplon with wave 
number q. The momentum is transferred to the electron, 
and no transitions between electron states occur. We will 
consider terms with q 7^ q'; the term with q = q' gives 
the shift of the electron energy levels. 

The matrix elements Wqq'j are linear in the param- 
eters of the direct two-ripplon coupling (|l^) and 
quadratic in the parameters of the one-ripplon coupling 



(j,0,0|F_q,q-|j,0,0)e 



-q' -q ; 



(27) 



where Vq^™" 



q 



and 



q 



(j, 0, 0|K_qe~*''""|j', V, rriy). In calculating renormalization 
of the parameters Wqq'j due to one-ripplon coupling we 
disregarded the contribution from virtual transitions into 
different states of out-of-plane motion because they 
involve a large energy change (it is straightforward to in- 
corporate the corresponding terms). We also disregarded 
ripplon energies fiiOq compared to "hio^^ . 

A calculation similar to that for crystal lattice defects 



i gives the phase relaxation rate for the coupling 
3) in the form 



q.q' 

Xfi(a;q)[n(Wq') -|- l](5(LJq - Uqi), 



(28) 

where n(ijj) = [exp(?icj/fcBT) — 1]^^ is the Planck num- 
ber. It follows from Eq. (28) that only thermally excited 
ripplons with LOq ^ ksT /fi contribute to the rate F^. In 
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what follows we will estimate contributions to from 
different mechanisms of electron-ripplon coupling taken 
separately and will again ignore cross-terms, which con- 
tain products of coupling constants for different mecha- 
nisms. 

The contribution to the phase relaxation rate F^**^' from 
the direct two-ripplon kinematic coupling ( jl^ ) has a sim- 
ple form in the case where the frequencies of ripplons 
with q < l/a|| are small compared to ksT/h. Then 

F^'^'' is determined primarily by forward scattering of rip- 
plons off the electron, with |q — q'| ^ l/ayj but with 
LOq = ujg' ~ kgT/h 3> ujr- Calculating the integral 
over the angle between q and q' by the steepest descent 
method, we obtain 



J- A — 



27V2a„ V 



(29) 



where K12 is the difference of the expectation values of 
the kinetic energy pl/2m in the states 1 and 2 divided 
by R; we have K12 = 3/4 for E^^ — 0, and K12 decreases 
with increasing E±_ . The numerical value of F?^"* is very 



small for low temperatures, F 



(k) < 



0.7x10 



2 0-1 



forT=10 



mK and W||/27r = 20 GHz. 

The contribution from the direct two-ripplon polar- 
ization coupling ( p^ ) can be estimated by utilizing the 
fact that the wave vectors of thermal ripplons qt — 
{p/ay/^{kBT/h)^^^ are less than l/rs for low temper- 
atures. To lowest order in gr^s the polarization contri- 
bution is again given by Eq. (29), but now K12 is the 
difference of the expectation values of the potential en- 
ergy A/2z divided by R. The corresponding rate is of the 
same order as F^'^'' . 

We now estimate the phase relaxation rate due to one- 
ripplon coupling (|l0|). We note first that the kinematic 



terms in (10) drop out of the matrix elements Wqq'j, be- 



cause they do not have diagonal matrix elements on the 
functions | j, 0, 0) (j = 1,2). The terms quadratic in the 
electric field E±_ drop out from the difference Wqq'i— 'i'qq'2, 
because they are independent of the electron state normal 
to the surface. The major contribution therefore comes 
from the polarization one-ripplon coupling cx Upoi in 

(|l^). We will denote it as T^^°^\ 

Polarization terms in the operators do not depend 
on the in-plane electron coordinate. This makes it possi- 
ble to calculate the sum over z^, in Eq. ( p7| ) for Wqq' j . 
The calculation is simplified in the case ksT ^ cOr where 
q 3> l/a|| and q • q' ~ qq' S> a~^. Then the sum of 

VqqT^/Tii^Wii can be approximated by 
2(jlVLq|j)(j|V^q'b->exph(q- q')'a;;/4] {nu:M)-\ 



In this approximation we obtain 

'keT 



r(pol) P 



Tia 



R^k- 



12- 



(30) 



Here we assumed that the coefficient ki2 = 
|(l|z;poi((?z)|l)|2 - \{2\vpoiiqz)\2)\'^ is a smooth func- 
tion of q for actual q ^ qx- Its numerical value is 
^ 0.23 for q Ki qT and T=10 mK, it weakly depends 
on the pressing field Ej_. The phase relaxation rate is 

^0^°'^ ^ 10^ s^^- The one-ripplon polarization coupling 
is therefore a major mechanism of phase relaxation for a 
confined electron. 

The overall ripplon-induced phase relaxation rate ap- 
pears to be small. It displays an unusual temperature 
dependence T^, and comparatively weakly depends on 
the in-plane frequency . We note that it is much less 
than our previous estimate [ p^ obtained for the case of 
in-plane confinement by a magnetic field. 



RIPPLON-INDUCED SIDEBAND 
ABSORPTION 



An important consequence of coupling to ripplons is 
the occurrence of sidebands in the spectrum of microwave 
absorption by a confined electron. Sidebands are formed 
by an electron 1 1,0,0) |2,0,0) transition accompa- 
nied by creation or annihilation of one or several rip- 
plons. Ripplon sidebands for a confined electron are sim- 
ilar to phonon sidebands in absorption spectra of defects 
in solids [Q. They can be understood from the Franck- 
Condon picture of an electron transition as happening 
for an instantaneous ripplon configuration. Since equi- 
librium ripplon positions are different in the ground and 
excited electron states, the transition is accompanied by 
excitation or absorption of ripplons, and the transition 
energy differs from its value E2 — Ei in the absence of 
coupling to ripplons. 

In order to describe the effect it suffices to keep only 
diagonal in the relevant electron states \j, 0, 0) part of 
the Hamiltonian of electron-ripplon coupling. One can 
then apply a standard canonical transformation which 
shifts ripplon coordinates so that they are counted off 
from their equilibrium values in the ground electron 
state. The transformed one-ripplon interaction Hamil- 
tonian (^, (^^ then takes a Franck-Condon form 

- E eqA^^('z)|2,0,0)(2,0,0|, (31) 
F{q) = <z2[^2|z;poi(9^)|2)-(l|«poi(9^)|l)]e"''""^'. 

For weak coupling, of primary interest are one-ripplon 
sidebands. Because ripplon occupation numbers are big 
for fc^T 3> hujr, the absorption cross-sections for electron 
transitions accompanied by absorption and emission of a 
ripplon are the same. Respectively, the sidebands are 
symmetrical as functions of frequency detuning Auj = 
Lj-{E2-Ei)/h (wehave|Aw| ~ < {E2-Ei)/n). The 
absorption is quadratic in the electron-ripplon coupling 
parameters and can be calculated by perturbation theory 
in Hf''^. From (|3l|) we obtain for the scaled sideband 
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FIG. 5: Scaled absorption coefficient in the sideband ( [3^ ) 
vs. frequency detuning Au — uj — {E2 — Ei)/h for E± = 
and uj'/2n = 20 GHz. 



absorption coefficient ashi^) 

■kTi ct^Oii 

a,b = a,| j dqq-^F^{q)6{Au; ± tUq). (32) 

The full absorption coefficient is given by ash multiplied 
by the integral over frequency of the zero-ripplon absorp- 
tion coefficient. The halfwidth of the zero-ripplon line T 
is the sum of the decay and dephasing rates calculated 
above. In ( ^2|) we assumed that \Auj\ ^ T; this inequal- 
ity is well satisfied in the interesting region \Auj\ ^ Uj., 
since from the above estimates a^^/F > 10"*. 

The intensity of the sideband is determined by the fac- 
tor Gsb- For T=10 mK and lj^^/2tt = 20 GHz we have 
Gsh ~ 0.1. The smallness of Gsb indicates that side- 
bands formed by two- or many-ripplon processes are not 
important. 

The scaled absorption coefficient in the one-ripplon 
sideband is shown in Fig. ^. It monotonically de- 
creases with the increasing distance |Aw| from the zero- 
ripplon line. For small \Auj\ (but \Auj\ ^ F) we 
have ash oc |Aw|~^/'^. As expected, decay of the side- 
band with increasing \Alu\ is much slower than de- 
cay of the Lorentzian tail of the zero-ripplon line oc 
r/(Aw)^. For large |Acj|/cJr, the sideband absorption 
falls off as exp[—{\ALu\/LUr)'^^'^ /2], because coupling to 
short-wavelength ripplons is exponentially weak. We 
note that the one-ripplon sidebands do not display struc- 
ture, in contrast to sidebands in electron-phonon systems 
in solids that reflect singularities in the phonon density 
of states. 



A. Intensity of the zero-ripplon line 

The integral intensity of the electron absorption spec- 
trum (integral over frequency of the absorption coeffi- 
cient) is determined by the matrix element of the dipo- 
lar transition cx Z12 and is independent of the electron- 



1.2 
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eEn,/R 

FIG. 6: Scaling factor W in the expression for the Debye- 
Waller exponent W. The curves 1 to 4 refer to rs/ay = 
0.18, 0.25, 0.35, 0.5 (l^,, /2-k ^ 10, 20, 39, 79 GHz, respectively). 

ripplon coupling for E2 — Ei 3> hujj.. However, the in- 
tensity of the zero-ripplon line is reduced, because of the 
sidebands. This reduction is described by a Debye- Waller 
type factor (or Pekar-Huang- Reese factor in the theory of 
electron-phonon spectra) exp(— FF). The parameter W 
is given by the integral of ctsb over to, 

W ^GshW,W = 2a^^ j dqq-^F^(q). (33) 

The dependence of the scaling factor W on the field 
E±_ and a^is shown in Fig. ^. It is clear from this figure 
and Eq. (B3|) that W weakly depends on the in-plane elec- 
tron frequency lo^^ as long as the corresponding ripplon 
frequency <C kBT/fi. At the same time, W decreases 
with the increasing pressing field E± , because the differ- 
ence in the effective radii of the electron states |1) and 
1 2) decreases, and so does the difference in the ripplon 
equilibrium positions in the states |1) and |2). 

The calculated value of the Debye- Waller factor for 
electrostatically confined electrons W w 0.1 — 0.05 for 

/27r = 20 GHz) and E^ varying from to - 300 V/cm) 
is close to the estimate W ~ 0.05 given earlier ||l^ for 
the case of in-plane confinement by a magnetic field. This 
factor emerges also in the analysis of the operation of a 
quantum computer based on trapped atomic ions [ p8[ , 
because optical transitions are connected to vibrational 
modes of the ions (the number of such modes is small, 
for a small number of ions) . 

In the context of quantum computing, sideband ab- 
sorption and the Debye- Waller reduction of the zero- 
ripplon absorption strength differ qualitatively from elec- 
tron decay and dephasing. The Debye- Waller mechanism 
does not affect an electron qubit between quantum op- 
erations. In contrast to dissipative effects, it does not 
accumulate between operations. However, it shows that 
a fraction of electron transitions may go wrong, as they 
are accompanied by excitation of ripplons. The num- 
ber of such transitions, and therefore the role of the 
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Debye- Waller factor, can be significantly reduced by us- 
ing longer excitation pulses, as we will discuss in a sepa- 
rate publication. 

VI. DECAY AND DEPHASING FROM 
COUPLING TO THE ELECTRODE 

Coupling to the underlying electrode may also provide 
an important mechanism of relaxation of a confined elec- 
tron (qubit). The corresponding relaxation parameters 
can be analyzed in a standard way. Fluctuations of the 
electrode potential modulate the inter-level distance and 
thus give rise to dephasing. In addition, an electron can 
make a transition between the states, with energy being 
transferred to an excitation in the electrode (for exam- 
ple, an electron- hole pair). The analysis is simplified by 
the fact that the size of the wave function of the qubit 
^ tb is small compared to the distance to the electrode h. 
Then the interaction with the electrode can be described 
in dipolar approximation, 

Hdip^~eS£±z, (34) 

where S£± is the fluctuating part of the field on the elec- 
tron normal to helium surface. This field comes from 
charge density fluctuations in the electrode. Here and 
below we do not consider effects of fluctuations of the 
electrode potential on in-plane motion, which are weak 
and not important for qubit dynamics. 

Electron relaxation parameters can be expressed in 
terms of the correlation function of the fluctuating fleld 

/•OO 

Q{oj)= dte"^\S£±{t)S£±{0)). (35) 
Jo 

We will assume that the function Q{u!) is smooth in the 
frequency ranges of interest for dissipation effects, i.e. for 
either tu < ksT/h or w ~ (£'2 - Ex)/n. 

Field-induced time variation of the phase difference of 
the wave functions |1) and |2) comes from Stark shift and 
is equal to 

Classical (thermal or quasi-thermal) field fluctuations 
give rise to phase diffusion on times that largely exceed 
the correlation time of S£±{t), which we assume to be 
short, < h/ksT. The coefficient of phase diffusion is 

equal to the dephasing rate F^'^''' . From Eq. ( ^Sf ) 

T^;'^ ^e^{z22~ziif'^eQ{Q)/h\ (36) 

If the noise spectrum Re Q{'jj) has peaks at low frequen- 
cies w < r^°'\ or the noise 5£±{t) is non-Gaussian, decay 
of (exp[i(5(/)i2(i)]) becomes nonexponential. Although the 
analysis has to be modified in this case, it is still conve- 
nient to relate decoherence of electrons on helium to the 
fiuctuating field 5£±{t). 



Decay rate of the qubit is determined by the prob- 
ability of a field- induced transition |2) — |1) between the 
electron states. This probability is determined, in turn, 
by quantum fluctuations of the field 5£^{t) at frequency 
^12 = {E2 - Ei)/h. From Eq. (H), 

ris" = e2|zi2pRe Q{n^2)/Ti\ (37) 

Here we assumed that decay is due to spontaneous emis- 
sion only, i.e. that there are no induced processes with 
energy transfer E2 ~ Ei. 

To estimate relaxation parameters we will assume that 
the controlling electrode is a conducting sphere of a small 
radius rd submerged at depth h beneath helium surface, 
as discussed in Sec. II. For low frequencies the surface of 
the sphere is equipotential. Then the fiuctuating field of 
the electrode is simply related to its fluctuating poten- 
tial SVci, S£± = 5Vc\rc\/h?- Much of low-frequency noise 
is due to voltage fluctuations from an external lead at- 
tached to the electrode, which has resistance TZcxt and the 
temperature Text that largely exceeds the helium temper- 
ature T. The noise can be found from Nyquist's theorem 
and gives the dephasing rate 

Ff'^ - 2fcBText7^cxte'(z22 - z^ifrlJh^h\ (38) 

For 7?.cxt = 25 Ohm, Text = 1 K, ^ci = 0.1 /im, h — 
0.5 /xm, and Z22 — ^ii = rs we obtain F^°''' w 1 x 10"^ s~^. 
This shows that thermal electrode noise may be a major 
source of dephasing for a qubit. The requirement to keep 
this noise small may be important in determining the 
depth by which controlling electrodes can be submerged 
below helium surface. 

In contrast to low-frequency noise, high-frequency volt- 
age fluctuations from sources outside the thermostat can 
be filtered out. Much of high-frequency quantum fluc- 
tuations that affect a qubit come from the underlying 
microelectrode itself. They depend on the interrelation 
between the electron relaxation time Tei in the electrode 
and ^12 ■ If "^01^^12 ^ 1, the electrode conductivity does 
not display dispersion up to frequencies ^ ^12] it greatly 
exceeds Q.12 for typical Vti2. 

An order-of-magnitude estimate of the decay rate r^2''' 
can be made by assuming that the controlling electrode is 
a lead attached to a sphere, and this sphere is equipoten- 
tial (fiuctuations of the total charge in the sphere make 
a major contribution to the field 6£±_ for small rei/h). 
Then from Nyquist's theorem 

r(°') = 2{E2 - E,)n,ye'\z^2?rl/h^h\ (39) 

where TZc\ is the resistance of the lead. If we estimate 
it as 0.1 Ohm, then using the same parameters as in 
the estimate of F^*^'' and setting E2 — Ei equal to the 

"Rydberg" energy R @), we obtain F^2^ ~ 5 x 10^ s'^. 
Even though this estimate is very approximate, it is clear 
that the major effect of electrodes on qubit relaxation is 
dephasing rather than decay. 
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VII. CONCLUSIONS 

In this paper we have provided a quantitative anal- 
ysis of the parameters of qubits based on electrons on 
helium. We introduced a simple realistic model of elec- 
trodes, which are submerged into helium in order to lo- 
calize and control the electrons. This model allowed us to 
estimate parameters of the electron energy spectrum and 
their dependence on the electrode potential. Control is 
performed by varying the field E± normal to helium sur- 
face. This field changes the distance between the energy 
levels of a qubit, which are the ground and first excited 
levels of motion normal to the surface, and enables tuning 
qubits in resonance with each other and with externally 
applied microwave radiation. 

The electrode potential determines not only E±, but 
also the in-plane electron confinement. We found the fre- 
quency of electron vibrations parallel to helium sur- 
face and related it to the field E±. Typical frequencies 

/27r are of order of a few tens of GHz, and typical fields 
are ~ 100 — 300 V/cm. We analyzed both the cases of 
one electrode and an electrode array, and investigated 
the effects of electrode geometry, including the inter- 
electrode distance and the depth by which electrodes are 
submerged into helium. 

We identified relaxation mechanisms, estimated decay 
rates for a confined electron, and found their dependence 
on control parameters. In contrast to unconfincd elec- 
trons studied before, decay is due primarily to electron 
transitions in which energy is transferred to two ripplons 
propagating in opposite directions or to a bulk phonon 
propagating nearly normal to the surface. In both cases 
helium excitations with comparatively large wave num- 
bers are involved. For different coupling mechanisms we 
have found the dependence of the decay rate on the pa- 
rameters of a confined electron. The decay rate is essen- 
tially independent of temperature, for low temperatures. 

The overall decay rate is of order 10^ s~^ for typical cu^^ . 
This estimate is obtained assuming that the typical wave 
numbers of excitations into which an electron may scatter 
are ^ 10^ cm~^. We expect that coupling to ripplons and 
phonons with much shorter wave lengths is small. Then 
the decay rate can be significantly decreased if a magnetic 
field is applied perpendicular to helium surface, because 
such field leads to an increase in level spacing of in-plane 
electron excitations, and therefore more energetic helium 
excitations have to be involved in decay. 

The major mechanism of dephasing due to coupling 
to excitations in helium is scattering of thermal ripplons 
off an electron. We calculated the scattering rate and 
showed that it displays an unusual temperature de- 
pendence. The major contribution to the dephasing rate 
comes from processes which involve virtual transitions 
between electron states. The ripplon-induced dephasing 
rate is ^ 10^ s~^ for typical lu^^ and T = 10 mK. 

An important mechanism of dephasing is voltage fluc- 
tuations of controlling electrodes. The dephasing rate 
strongly depends on the source of these fluctuations and 



also on the depth by which electrodes are submerged 
into helium. An estimate for Johnson noise from a typi- 
cal lead connnected to an electrode gives dephasing rate 
- 10^ s-\ 

We have also analyzed sidebands of electron absorp- 
tion spectrum related to electron transitions accompa- 
nied by emission or absorption of a ripplon. We found the 
Debye- Waller factor which describes the intensity of the 
zero-ripplon absorption line and characterizes the over- 
all probability of exciting a ripplon during an electron 
transition. This factor gives fldelity of qubit excitation 
by microwave radiation, which is a major single-gate op- 
eration of the quantum computer based on electrons on 
helium. 

The results provide a quantitative basis for using elec- 
trons on helium as qubits of a quantum computer. The 
clock frequency of such computer flqc, which is deter- 
mined by the dipole-dipole inter-electron interaction, is in 
the range of 10^ — lO'* Hz even for inter-electron distances 
« 1 /im, and therefore it largely exceeds both decay and 
dephasing rates of a confined electron. Our results sug- 
gest ways of further reduction of these rates. They show 
how to choose parameters of the system in an optimal 
way and also show that there is an extremely broad range 
where the parameters can be dynamically controlled, be- 
cause the inter-level distance E2 — E1 huj^^ fiV,Qc ^ 
nr. They also suggest a sequence of steps that have to 
be done in order to implement a quantum computer with 
electrons on helium in experiment. 



Acknowledgments 

We are grateful to A. Dahm, B. Golding, and J. Good- 
kind for valuable discussions. This research was sup- 
ported by the NSF through grant No. ITR-0085922. 



APPENDIX: ONE-RIPPLON POLARONIC 
EFFECT 

Besides relaxation, coupling to ripplons leads also to 
a polaronic effect. Because ripplon frequencies are low, 
the major contribution comes from processes in which 
a ripplon is created or annihilated, but the state of the 
electron system is not changed. Polaronic shift of the 
electron transition frequency is then determined by the 
diagonal matrix elements of H^^^ on the wave func- 
tions |1,0, 0), |2,0, 0). Keeping only these terms in H-^^ 
corresponds to the adiabatic approximation in which rip- 
plons have different equilibrium positions depending on 
the presence of an electron (one can think of a "dimple" 
made by an electron on helium surface pj|) and on the 
electron state. Of primary interest to us is the state de- 
pendence, as it characterizes the strength of coupling of 
the electron transition to ripplons. The corresponding 
coupling is described by the Franck-Condon interaction 
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FIG. 7: The factor fp in the Franck-Co ndon polaronic shift 
of the transition frequency of a qubit (A.l) as a function 
of the pressing field E± for typical values of the in-plane 
localization length a^^. The curves 1 to 4 correspond to 
tb/o,^^ = 0.18,0.25,0.35,0.5; the respective values of uj^J2-k 
are ^ 10, 20, 39, and 79 GHz. 



Hamiltonian (13]]) 



The Franck-Condon polaronic shift of the transition 
frequency 1 — > 2 due to the coupling ( |3l| ) is given by a 
simple perturbation theory, 



P = Pi2fp{a^^,E^_), Pi2^RV4Arl, (A.l) 
Jo 



where fp is a dimensionless factor determined by the ma- 
trix elements of Vpoi (10) on the wave functions of out-of- 
plane motion. It depends on the dimensionless parame- 
ters ol^Jtb and eExTsIRi and is numerically small for 
typical parameter values, see Fig. |^. 

The numerical value of the factor P\2 is Pi2/27r « 
2.2 X 10^ Hz. The energy TiPyi is much less than the 
distance between the electron energy levels. The shift P 
is also smaller than the typical frequency of ripplons 
coupled to the electron (p5|). The inequality \P\ <C 
indicates that the |1) |2) transition is weakly coupled 
to ripplons. 
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